Quantum-fluctuation effects on the thermopower of a single-electron transistor 
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We study thermal conductance and thermopower of a metallic single-electron transistor beyond the limit of 
weak tunnel coupling. Employing both a systematic second-order perturbation expansion and a nonperturbative 
approximation scheme, we find, in addition to sequential and cotunneling contributions, terms that are associated 
with the renormalization of system parameters due to quantum fluctuations. The latter can be identified by their 
logarithmic temperature dependence that is typical for many-channel Kondo correlations. In particular, the 
temperature dependence of thermopower, which provides a direct measure of the average energy of transported 
particles, reflects the logarithmic reduction of the Coulomb-blockade gap due to quantum fluctuations. 

PACS numbers: 73.23.Hk, 73.50.Lw, 85.80.Fi 
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I. INTRODUCTION 

Transport of electrons through a small metallic island is 
strongly affected by charging effectsi 2 ^ Tunneling of an 
electron on an island with capacity C is associated with an en- 
ergy of the order of the charging energy Ec = e 2 /2C. At low 
temperature, kgT <C A, where A is the charging-energy gap 
between ground state and first excited charge state, transport 
is suppressed. In a single-electron transistor (SET), an island 
connected to two leads by tunneling junctions (see Fig.^, this 
blockade of transport can be controlled by an additional gate, 
resulting in the well-known Coulomb oscillations of current 
with respect to gate voltage4^ 

If the island is well isolated from the leads, i.e., the barrier 
resistances R^ R are high, 



fects such as the thermal conductance and the thermopower 
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electric transport is dominated by first-order transport in 
the tunnel conductance ao (sequential tunneling). In the 
Coulomb-blockade regime, where sequential tunneling is ex- 
ponentially suppressed, inelastic cotunneling becomes im- 
portant. In these processes of second order in ocq, the en- 
ergetically unfavorable charging of the island occurs only 
virtuallypli 7 ^^ But also at resonance, where sequential tun- 
neling is present, there are higher-order transport contribu- 
tions. They are associated with renormalization of charg- 
ing energy and tunnel conductance due to quantum fluctua- 
tions. This can be qualitatively understood by mapping the 
SET at low temperature and close to resonance to a many- 
channel Kondo problem and performing a poor man's scal- 
ing analysis of the latter. 10 For a quantitative analysis of these 
quantum-fluctuation effects, a systematic second-order pertur- 
bation expansion within a diagrammatic real-time technique 
has been performed 1112 and used to study different single- 
electron systemS(A2iA2iA£ i n particular, a logarithmic reduc- 
tion ~ Ooln fiEc of the maximum conductance, indicating a 
renormalization of the tunnel conductance, has been found, in 
quantitative agreement with experimental observationsiliiii 

In electron transport, the transfer of charge and heat are 
connected to each other. This gives rise to thermoelectric ef- 
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where V is the voltage due to a temperature difference 8T in 
the absence of a charge current /. With the help of Onsager 
relations, 18 the thermopower can be related to the average en- 
ergy (e) of the transported electrons relative to the Fermi en- 
ergy: 



S=- 



eT 



(2) 



In macroscopic conductors, the thermopower is of the order 
of {kg / e)(kBT / Ep), which, in general, is very small. This is 
due to the fact that the product of density of states and elec- 
tron velocity squared, determining the contribution of elec- 
trons of a certain energy to the current, varies only slowly, 
namely on the scale provided by the Fermi energy. There- 
fore, one can perform a Sommerfeld expansion of (e). The 
zeroth-order term of the Sommerfeld expansion vanishes. The 
next-order correction, that accounts for an asymmetry of the 
product mentioned above around the Fermi energy, yields an 
extra factor of kgT /Ep. This is different in mesoscopic sys- 
tems for which charging effects can strongly modify the effec- 
tive density of states. In the cotunneling regime of a SET, for 
example, the effective density of states changes on the scale 
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FIG. 1: Setup for thermopower measurement on a single-electron 
transistor. The two leads are kept at different temperatures and a 
voltage bias V = Vl — Vr is applied, such that no electrical current 
flows through the device. A gate voltage sets the working point and 
thereby the charging-energy gap A of the SET. 
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A, i.e., the thermopower is of the order of (kg / e) (IcbT / A) . An 
even larger thermopower is generated by sequential tunneling 
processes. They are modeled by a delta-function peaked den- 
sity of states around A, for which the Sommerfeld expansion 
is not applicable. Instead, a direct evaluation of (e) ~ A leads 
to S ~ (fcfl/e)(A/fcfl!T). In conclusion, at the crossover from 
sequential to cotunneling the thermopower can reach values 
of the order of ks/e. 

By now a few experiments on thermal conduc- 
tance or thermopower in quantum dots have been 
performed. 19 - 20 - 21 ' 22 - 2 i2i2^ Thermoelectric effects in 
various mesoscopic systems have also been studied 
theoreticallyi 26 i 27 i 28 i 29 i 3Q i 3 1 1 32-,33 1 34 1 35 Q uantum dots with 
discrete 20 (single-particle) energy level spectrum as well as 
dots closer to the metallic (quasi-continuous) limit 19 have 
been compared to sequential-tunneling theory. For increased 
coupling and lower temperatures, systems of discrete level 
structure showed signatures of cotunneling 21 and also of 
Kondo physicsi 2 ^ For the system under consideration in 
the present paper, a metallic single-electron transistor, 
thermopower has been investigated by taking into account 
sequential 26 and cotunneling 27 processes. 

In this paper, we study the effect of quantum fluctuations 
on the thermal conductance and thermopower of a metallic 
single-electron transistor with a large number of transverse 
channels. Thermopower provides a direct access to measure 
a renormalization of the charging-energy gap due to quantum 
fluctuations since it is closely related to the average energy of 
tunneled electrons. As we predict below, this will give rise to a 
logarithmic temperature dependence of both the slope of ther- 
mopower at resonance and of the position of the maximum of 
thermopower as a function of gate voltage. 

This work is structured as follows. First, in Sec. [[I] we de- 
fine the model Hamiltonian of the system under consideration 
and derive expressions for the thermal conductance making 
use of two different approximation schemes. In Sec. Mil we 
discuss how the thermal conductance and the thermopower 
exhibit the renormalization of the charging-energy gap in- 
duced by quantum-fluctuations. We summarize our results in 
Sec.HVl 



II. THEORY 

A. The system 

A metallic single-electron transistor (see Fig.Q is modeled 
by the Hamiltonian 

H = H L +H R + H, + H ch + H T =H + H T . (3) 

Here H r = Zk„ 4n a lkn a rkn and Hj = Y. qn Zq n c\ n c qn describe 
noninteracting electrons in the two leads r — L,R, and on the 
island, respectively. The index n = 1 , . . . , Nt is the transverse 
channel index which includes the spin while the wave vec- 
tors k and q numerate the states of the electrons within one 
channel. In the following, we assume the many-channel limit 
ffi>l. Coulomb interaction of the electrons on the island is 



described by the capacitance model =Ec(n — n x ) 2 , where 
Eq = e 2 1 (2C) with total island capacitance C = Cl + Cr + C h . 
This electrostatic energy depends on the number of excess 
electrons on the island, given by their number operator h, as 
well as on applied gate and bias voltages. The latter are ac- 
counted for by the "external charge" en x = ClVl + CrVr + 
CgVg. To increase the number of electrons in the island 
from N to N + I one has to overcome the charging-energy 
gap, the difference between neighboring charge states A^ = 
(N+l\H ch \N+l) - (N\H ch \N) =E C [1 + 2 (N-n x )}, which 
is tunable by the gate voltage via n x . The resonance condition 
An = 0, where the charging-energy gap vanishes, is fulfilled 
at half-interger values of n x . 

Finally, charge transfer processes are described by the tun- 
neling Hamiltonian 

H ?= I W4knC< ! ne- i *+h.C.. (4) 
r=L,Rkqn 

The matrix elements Tj™ = T r are assumed to be independent 
of the states k and q and channel index n. They determine 
the tunneling resistance Rj, r of the left and right junction by 
l/R T , r = {2ne 2 /h)N t N r (0)N I (0)\T r \ 2 , where N I/r (0) are the 
density of states of the island/leads at the Fermi level. Note 
that, while the number A, of channels is large, the contribution 
of any one channel is so small that the total coupling remains 
weak and Coulomb blockade effects will occur. The operator 
e ±l(? shifts the charge on the island by ±e. Since the left and 
right lead are assumed to be reservoirs with fixed electrochem- 
ical potential ji r and temperature T r , the phase <p (or its canon- 
ical conjugate, the island charge n) is the only independent 
dynamic variable in our model. In general, the electron tem- 
peratures of left lead, island, and right lead can all be different 
from each other and differ from the lattice temperature. 36 

B. Conductance, thermal conductance, and thermopower 

The current I r flowing into reservoir r can be expressed by 
using correlation functions for the island charge C > (f,f') — 
-i(e-m) e m')) and C<(f/) = j^'V^')). For a time- 
translational invariant system these correlation functions de- 
pend only on the time difference, C(t,t') — C(t — t'), and we 
will work with the Fourier transforms C(co) = J dte tmt l h C{t). 
The tunneling current / = //, = —Ir is determined by 

ie f 

I r = -- dco[a r+ {(o)C > {co) + a r -{co)C < {co)} , (5) 

which includes all possible tunneling processes via the exact 
correlation functions C<(o)). The rate functions 

a r± {co) = as£dEf?(E + (o)f*(E) (6) 

enter, a r+ (co) describing tunneling of an electron from lead r 
onto the island, and a r ~(co) from the island to lead r. Here, 
/+ denotes the Fermi function, and /~ = 1 — f + . Applied 
temperature or voltage gradients, 8T = Tl — Tr and V = Vl — 
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Vr, are accounted for by evaluating f^(E + o) at temperature 
T r = T + 8T, and voltage V r , while f^(E) is taken at T. We 
define also a(to) = £ r a , ' + (a)) + a r ~(a)). 

The linear electrical and thermal conductances are given by 



dV 



and 



Gr = 



37 



V,ST=Q 



dST 



(7) 



v,5r=o 



respectively. The thermopower describes the voltage gener- 
ated by a temperature difference in the absence of an electrical 
current, and is related to the above mentioned conductances 
by 



lim -=— 

ST->0 OT 



7=0 



Gj 

g7 



To calculate the linear electrical and thermal conductance 
Gv and Gt, we expand the rate functions up to linear order 
in either V r or 8T r in the expression for the current, Eq. (|5jl- 
It is convenient to use current conservation £ r I, = to write 
the current as the combination I = (oCqIl — OCqIr) /(Oq + Og ). 
When expanding this combination up to linear order in either 
V r or 8T r , we immediately see that only the equilibrium cor- 
relation functions C^(to), taken at ST = and V — 0, enter, 
since linear corrections in V r or 8T r drop out in the combi- 
nation considered. In equilibrium, the correlation functions 
are related to the spectral density A(to) for charge excita- 
tions on the island by C^to) = -2ni[l — /(<o)]A((») and 
C<(co) = 2nif(to)A{to). This, eventually, leads to the linear 
electrical and thermal conductance 



and 



Gv — G a 



Gt — —Gas — 



sinhj3<a 



dtO^A(to), 



sinh/3a> 



(8) 



(9) 



respectively, where G as = l/(Rj+Rf) is the classical 
electrical conductance asymptotically reached in the high- 
temperature limit. In conclusion, we need to evaluate the 
equilibrium spectral density A((o) to obtain the linear elec- 
trical and thermal conductance via Eqs. ^ and l|9}. To keep 
notation simple, it is convenient to introduce dimensionless 
conductances gv and gr, defined by 



gv 



Gv , e Gt 

— and gr = -7-7^- (10) 

Gas KB Gas 

As we see from Eqs. (|SJl and the dimensionless conduc- 
tances differ from each other by a factor j3 to/ 2 in the inte- 
grand. This can be easily understood with the help of Eq. (|2ji, 
which indicates that to/2 is the average energy of the lead 
electrons (measured relative to the Fermi energy) that con- 
tributes to a island charge excitation of the energy to. The fac- 
tor 1 /2 comes from averaging over the available phase space 
of the electronic states in the leads and the island. Roughly 
speaking, on average one half of charge excitation energy 
comes from or goes to the lead and the island electrons, re- 
spectively. 



C. Approximation schemes 

In the following we will employ two approximation 
schemes for calculating the spectral density and, thus, the 
linear electrical and thermal conductance. On the one hand, 
we will perform a systematic perturbation expansion up to 
second order in the dimensionless tunnel conductance ao- 
On the other hand, we will use a nonperturbative resum- 
mation scheme, the so-called "resonant-tunneling approxi- 
mation" discussed further below. Both schemes go beyond 
the weak-coupling (sequential-tunneling) limit of small tun- 
nel conductances, but in different ways. The virtue of either 
scheme as compared to the other one is discussed below. Both 
of these schemes are based on a real-time diagrammatic tech- 
nique introduced in Ref.[^[ Here, we will make use of known 
results of these methods without the need for an explicit recal- 
culation of the diagrams. Therefore, in this paper, we will not 
discuss rules for constructing and evaluating diagrams, but re- 
fer the interested reader to the existing literature. 



1. Systematic perturbation expansion 

We perform a systematic perturbation expansion of the 
correlation functions C^(to) = {®) an d> there- 

fore, automatically for the spectral density A(to) = [C < (to) — 
C > (to)]/(2ni) = Y£ = oA( k \tQ), where the index k denotes the 
power of Oo in the expansion. The real-time method yields 
diagrammatic representations of the correlation functions in 
different order (see Fig. 3 of Ref. fl2h . which are calculated as 
sketched in Sec. Ill of Ref.E 

To lowest order, the spectral density needed is simply 



A se< !(a)=£(/V 

N 



-P N+ i)8(to-A N ) 



(11) 



with the equilibrium probabilities (to zeroth-order in c«o) 
P N = exp[-j5E ch {N)}/Z with Z = £exp [-j5E ch (N)} , 

N 

to find the island in charge state N. As a result, the dimension- 
less linear electrical and thermal conductances are 



seq 
#7 



N 



N 



nv+i 



sinh/3A;v 

(^A w /2) 2 
sinhj3Aft( 



(12) 
(13) 



For low temperatures at most two charge states contribute, 
e.g., for n x w 0.5 only the term YV = enters, and Pq+P\ = I. 
Since, in the lowest order, the only allowed charge excitation 
energies are Ajv, the average energy of the contributing elec- 
trons for transitions between charge state iV and iV+ 1 is An/2. 

For the next-order contribution, we use correlation func- 
tions in the limit of vanishing applied voltage and temperature 
bias from Ref. [3 namely Eqs. (12), (14), and (15) together 
with Eqs. (20) and (60) there. We can perform all integrals in 
Eq. (|9j above analytically to find the complete second-order 
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contribution to the linear electrical and thermal conductances 
as a sum of four terms, 

S?J T =^v t /T + 8v/T+^/T + 8v/T ■ (W) 



For the dimensionless thermal conductance, we get the 
lengthy but complete expressions 



cot 

8t 



g? 



gf 



8 2 r e 



N 



2 2 dN + a N -i <j>N-(j) N -i+A N -\d(j) N -\-A N d(j) N 
a N -\A N -id (p N -i+a N A N d <j) N -\ 



E c 



*tf sinhj3AA, 



d (2<p N + 0jv_i + <j> N+ i) + 



<j>N-l - 0N+1 



(14a) 
(14b) 



N 



j3E^V'(0iV'-i-0A") 
A" 

PA N /2 



PPn (<Pn-i - <Pn) + PPn+i - 



sinh/3A/\ 



Pn+Pn+i 

(P N + P N+ l) (20at - 0jv_i - 0jv+l) 



La N + a N -\ J3(A W + A w _i) ~~ \ 
o • ufl/A — T^v V v^- 1 +/ w+i) 
w 2 sinhj3(A iV + A w _i) 



2£ c 



Aiv+AAr_i 



(14c) 
(14d) 



r 



Here, <9 stands for d/dA^ = —(1 /2Ec)(d / dn x ) and we used 
the definition <p N — (Oq + 0!g )A w Re l P (ij3A N /2n), where *F 
denotes the digamma function. Furthermore, we defined 
aw = J3Ajv/2. The result for the dimensionless electrical con- 
ductance gv is the same but without the factors a^, a^— l an d 
(oat + aa?_i)/2 in accordance with Ref. 12, The factors 
and ajv_i account for the average energy of the lead electrons 
contributing to transport. 

For the interpretation of the four terms we follow the rea- 
soning put forward in Ref. 12, The first term gf l models co- 
tunneling processes, where an electron is transferred through 
the whole device without changing the charge of the island. 
This is the dominant transport contribution far from the reso- 
nance of sequential tunneling. We can identify this term with 
the "reguralized" cotunneling result postulated in Ref. |27j In 
fact Eq. ( I Hat stems from the integral expression 



N 



dCO CO 



(JSco/2) 2 



xRe 



2sinh 2 ()3co/2) 
1 



fi5ah 



co - A N + i'O 4 



CO 



-i0+ 



which yields Eq. (30) of Ref. 27 in the regime considered 
there. Note that the infinitesimal imaginary parts of the de- 
nominator arises naturally within the diagrammatic theory, 
not requiring regularization by hand as in Ref. |27|. At low 
temperature and away from resonance, e.g., in the Coulomb- 
blockade valley with Pq = 1, we can make use of the expan- 



sion Re l P(k) 



cot 

8v 



cot 

8t 



ao- 



= In \x 
2n 2 



l/(12x 2 
(k B T) 2 



ao — (k B T) 3 



1 

Ao 
1 

Ao 



1/(120/) 

A-i 
1 

A~i 



to get 



(15) 



— + — (16) 
Ao A_i ' 



in accordance with Ref. 1271 

Away from resonance, Sy'/j I s me on ly second-order con- 
tribution. When approaching the resonance, two more terms, 
Sy/j an d gy/T come pl a Y- They are associated with 
sequential-tunneling processes but with renormalized param- 
eters: gyj T is the first correction term to sequential tunnel- 
ing due to renormalization of the tunnel-coupling strength, 
gyj T the respective correction due to a renormalized charging- 
energy gap. The relation of these terms to renormalization is 
discussed in more detail below. At low temperature, only the 
term N = contributes and Pq +P\ = 1, so that Eqs. J14bl > and 
J14ci reduce to 



a _ Q3Ao/2) 

ST 



8t 



A =d 



sinhj3Ao 

(J3Aq/2) 



sinhj3Ao 



<9(20 o + 0_i + 0iH 
(20t,-0_i-0i) 



0-i-0i 
Ec 



Obi) 



The fourth term, gj e describes cotunneling processes in 
which the charge of the island is changed by 2e. Since the total 
change of the charging energy between charge state N+l and 
N — 1 is Ajv + Ajv_i, the factor (a# + ajv-i) /2 accounts for the 
average energy per contributing lead electron. To overcome 
the charging energy for two electrons, a large temperature is 
required. The term g^f vanishes at low temperature, and hence 
will not be of importance in the following. 

The virtue of the perturbation expansion lies in the fact 
that (i) all second-order contributions are systematically taken 
into account, (ii) their identification with cotunneling pro- 
cesses and renormalization corrections to sequential tunnel- 
ing is straightforward, and (iii) all expressions are unambigu- 
ously fixed by the system parameters without any remaining 
cutoff parameters. With increasing tunnel-coupling strength 
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or lowering temperature, however, the second-order perturba- 
tion theory will become insufficient. Therefore, we also apply 
a different approximation scheme as described in the the fol- 
lowing. 



2. Resonant-tunneling approximation 

The so-called resonant tunneling approximation (RTA) has 
been introduced in Ref. |37| as a nonperturbative treatment of 
quantum fluctuations. It amounts to resummation of a certain 
diagram class, including contributions of arbitrary high order 
in the tunnel-coupling strength. In particular, only two charge 
states N — 0,1, and only density-matrix elements that are at 
most twofold off-diagonal are taken into account. For details 
of the derivation, we refer to Ref. |^ 

Within RTA, the equilibrium spectral function is found^i to 

be 



with the self-energy 



a {03} 



a{co) = / dco' 



(q — Aq — a(co)\ 2 
a(co) 



(17) 



co - co' + i0+ 



Real and imaginary part are given by 

'PD 



Recr(<») = -2cxqC0 In 
Ima(co) = -na(co), 



2k 



where D is a high-energy cutoff of the order of the charging 
energy or band width. The expression for the electrical and 
thermal linear conductance follows from Eqs. (|8} and (|9). 

The virtue of the RTA is that, due to the resummation of 
higher-order contributions, lower temperature and higher val- 
ues of the tunnel-coupling strength are accessible. This is also 
indicated by the fact that the self-energy <j(co), describing 
renormalization of the charging energy-gap and the tunnel- 
coupling strength, appears in the denominator of Eq. d!7l >. 
On the other hand, the truncation of the Hilbert space to two 
charge states leaves a high-energy-cutoff dependence of the 
results. Thus RTA is suited for describing effects of the qual- 
itative temperature dependence due to quantum fluctuations 
at low temperature. For quantitative results at higher temper- 
ature, the systematic second-order perturbation expansion is 
more reliable. 



D. Renormalization effects 

The main result of this paper is the appearance of renor- 
malization effects in thermoelectric quantities. Therefore, we 
comment in this subsection on the relation between quantum- 
fluctuation induced renormalization and the electrical and 
thermal conductance in more detail. This discussion is equally 
valid for the results of the electrical and thermal conductance, 
and follows along the line of Ref.fT^ 



The notion of system-parameter renormalization is the cen- 
tral idea of all renormalization-group (RG) schemes. An ef- 
fective low-energy model is derived by successively integrat- 
ing out high-energy degrees of freedom in the leads. A poor 
man's scaling version of such an RG scheme for the two- 
charge-state approximation of the metallic SET has been per- 
formed in Ref. LUJ by mapping it to a many-channel Kondo 
model. During the RG procedure, both the tunnel coupling 
strength ccq and the charging-energy gap Ao becomes renor- 
malized as a consequence of the tunnel coupling between is- 
land and lead electrons. The renormalized values a and A are, 
within this poor man's scaling scheme, 



a 

CCq 



A 
Ao 



1 



l+2aoln(D/coc)' 



(18) 



where D is the high-energy cutoff (the smaller of charging 
energy or band width) and C0c the low-energy scale at which 
the RG procedure stops (here the larger of temperature IcbT or 
charging-energy gap Ao). As a consequence of the large num- 
ber of transverse channels in the tunnel contacts, the charging- 
energy gap and the tunnel-coupling strength are renormalized 
towards lower values, with a logarithmic dependence on the 
high- and low-energy cutoff. While the result is inherently 
nonperturbative (an expansion of the denominator yields all 
powers of ao), it is rather qualitative as the numerical coeffi- 
cient D/cOc is unknown. 

In the spirit of an RG picture, an effective low-energy the- 
ory of transport that takes into account renormalization is ob- 
tained by taking the sequential-tunneling formula but with 
renormalized system parameters a and A instead of the bare 
ones Oq and Aq. This amounts to 



G(ao,A ) = G seq (a,A) + regular terms , 



(19) 



for the electrical or thermal conductance. The "regular terms" 
represent higher-order contributions, such as cotunneling pro- 
cesses, that are not associated with renormalization. The latter 
are not included in the RG procedure, and are not considered 
in the following. 

To relate the second-order transport contributions to renor- 
malization, we expand Eq. ( 1191 up to second order in Oo, 



a 



G^(a,A) = — G s ^(ao,A ) 



(A-A ) 



dAn 



By comparison with Eqs. fll4bll > and dl4cll >. we obtain 



a 



do 



Ao 



1 -2oq< -1 +ln 



fpE C 



l-20o 



1+ln 



V n 



AoRe'P I i 



PAo 
2% 



Re*P i 



Mo 

2% 



Within the RTA we find21 



a 



A 
Ao" 



l+2oo 



In 



(20) 
(21) 



(22) 
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FIG. 2: Thermal conductance of an SET calculated to first and sec- 
ond order in tunneling coupling for ksT/Ec = 0.05 and Oq = 0.04. 
The inset shows the different contributions in second order, due to 
cotunneling (g™'), renormalization of the tunnel-coupling strength 

(gj) and of the charging-energy gap (gj). Renormalization leads to 
a broadening and a suppression of the thermal conductance as com- 
pared to the sequential-tunneling result and cotunneling yields an al- 
gebraically decaying contribution in the Coulomb-blockade valley. 



The RTA result is nonperturbative in Oo, and it resembles 
the structure of the poor man's scaling RG result in Eq. i ll M . 
Its numerical value remains undetermined as the exact form 
of the high-energy cutoff D is left unspecified. This contrasts 
to the result from second-order perturbation theory. There, 
all numerical constants are specified. On the other hand, the 
renormalization is determined only up to linear corrections in 
Of). Indeed, this correction can be considered as the lowest- 
order term of an expansion in 0C()ln(Ec / max{Ao,kgT}) [cf. 
Eq. (I22H . which serves as small parameter. We conclude with 
the remark that the interpretation of some of the second-order 
contributions as first-order transport but with renormalized pa- 
rameters was recently supported by analyzing the full count- 
ing statistics of electrical transport through a metallic SET. 38 
There, the functional dependence of the cumulant generating 
function on the counting fields enables an unambiguous iden- 
tification of sequential and cotunneling, in full support of our 
interpretation above. 



III. RESULTS 



A. Thermal conductance 

Figure [2] shows the first- and second-order contributions 
to the thermal conductance, i.e., the sequential-tunneling re- 
sult, Eq. ( I13> . and the different contributions to second order, 
Eq. dl4t . in the inset. For the low temperature considered here, 
gj e vanishes, as it corresponds to two electron processes — 
leaving the dot in a state, where the electrostatic energy is 
changed by more than Eq. As discussed above, gf l repro- 
duces the regularized cotunneling result; it is the dominating 



contribution away from resonance (j3Ao 3> 1), as it decays al- 
gebraically only. From Eqs. dl5l and d!6l >. we see that the 
cotunneling contribution to the electrical and thermal conduc- 
tance scales with (T/Aq) 2 and (T/Aq) 3 , respectively. The 
(T /Aq) 2 behavior of the electrical conductance is understood 
from the fact that each tunneling rate contributes a factor T 
while the A 2 , denominator is that of standard second-order per- 
turbation theory. At fist glance, one might expect a T /Aq be- 
havior for the thermal conductance due to the relative factors 
aw = PAfj/2 in Eq. d!4a> . However, it turns out that terms 
in the thermal conductance stemming from the lowest order 
in 1/ (j3Ao) in the expansion of O cancel out when expand- 
ing Eq. dl4a> . Since the expansion of *P(ix) has only even 
powers in x, the first nonvanishing contribution to the thermal 
conductance scales with (T/Ao) 3 . 

At resonance, the terms gj and gj, associated with renor- 
malization of the tunnel coupling a and the charging-energy 
gap A, become important. The renormalization of the tunnel 
coupling strength towards lower value results in a reduction of 
the peak height. The renormalization of the charging-energy 
gap shifts the system effectively closer to resonance and con- 
sequently yields a broadening of the resonance structure. In 
other words, the renormalization of coupling is reflected in the 
suppression of the maximum value of thermal conductance, 
the renormalization of the charging-energy gap in the shift of 
the maximum's position. 

The results for thermal conductance look rather familiar 
from the conductance results and do not clearly showcase un- 
expected features. Looking at the thermopower, however, we 
can gain new and interesting insights in the mechanisms of 
electron transport through our system and how it is influenced 
by quantum fluctuations. This is owed to the intuitively ap- 
pealing interpretation of thermopower as measure of the aver- 
age energy of transported particles, see Eq. @. 



B. Thermopower 

The thermopower as a function of the gate charge n x for 
different temperatures is displayed in Fig. [5] We show the full 
result (black lines) 



^seq 



,(2) 



Gseq 
v ■ 



(23) 



that takes into account all first- and second-order contributions 
to the electrical and thermal conductance, and, for compari- 
son, also the pure sequential-tunneling result S seq — G^ q / Gy q 
(gray lines). The thermopower vanishes at both integer and 
half-integer values of n x . At resonance, i.e., at half-integer 
values of n x , the thermal conductance vanishes due to a can- 
cellation of transport contributions from lead electrons above 
and below the Fermi level, that generate the same charge ex- 
citation [Fig. |3Ja)]. In the middle of the Coulomb-blockade 
valley, i.e., at integer values of n x , the zero is due to a can- 
cellation of processes that involve adjacent charge excitations 
[see Fig.|3jb)]. In between [situation sketched in Fig.|3Jc)], 
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FIG. 3: Thermopower within perturbative calculation for (Xq = 0.002 
and k B T/E c = 0.1 (solid line), 0.04 (dashed), and 0.01 (dotted). 



the thermal conductance and, thus, the thermopower is finite, 
with alternating sign at integer and half-integer values of n x . 

At larger temperature (solid line in Fig. [3}, sequential tun- 
neling dominates. Results for this regime have first been de- 
rived in Ref. |2fJ We recover a sawtooth behavior, with a lin- 
ear increase as long as transport predominantly involves only 
one transition N ^ N + 1 of charge states. Then the average 
energy of transported particles and correspondingly the ther- 
mopower is proportional to An/2. Around n x « N, the adja- 
cent transition iV — 1 <-» N, that contributes with an opposite 
sign, comes into play. This gives rise to a sharply falling edge 
of the sawtooth with a width given by temperature. 

At lower temperatures (dashed and dotted lines), sequential 
tunneling dominates transport only around half-integer values 
of n x , but cotunneling takes over in the Coulomb-blockade 
valley in between, resulting in a suppression of the rising edge 
of the sawtooth. Instead thermopower decays with T/A away 
from resonance, as seen from Eqs. (I15> and (I16> . 

These features of the thermopower have been explained by 
Turek and Matveev 27 by considering sequential plus cotun- 
neling processes (the terms g seq and g cot only). They postu- 
late a universal low-T behavior, whereby S scales as 5(J3 Aat). 
We find that this does not hold true for a complete higher- 
order calculation. Taking into account the renormalization 
of system parameters due to quantum fluctuations lets the 
thermopower deviate from universal behavior as illustrated in 
Fig. |4] in which the thermopower is plotted as a function of 
j3Ao. Conversely, these deviations allow an insight into the 
renormalization process and reveal the rich physics missed by 
taking into account cotunneling processes only. 

In the following, we will concentrate on two distinctive fea- 
tures of Fig. |4] On the behavior close to the resonance of se- 



quential tunneling (at half integer values of n x ) and on the 
position of the maximum of thermopower. 



1. Reduction of charging-energy gap 

Close to the resonance at n x = 1 /2, transport is associated 
with charge excitations 0^1, and thermopower is linear in 
n x . The sequential-tunnelig result S seq = — (&B/e)/3Ao/2 cor- 
responds to an average energy Aq/2 of the contributing lead 
electrons. Including cotunneling processes slightly reduces 
the slope (dotted line) by a factor that is independent of tem- 
perature. The slope is further reduced when taking into ac- 
count all second-order contributions. Using Eq. (|2jl as a def- 
inition of the average energy we see that the reduction of the 
slope reflects a reduction of the average energy of the con- 
tributing lead electrons. In Fig. [5] we display the average en- 
ergy defined via Eq. (0 close to resonance as function of tem- 
perature (this corresponds to the slope at n x = 1/2 in Fig.|4ji. 

Sequential tunneling gives the ratio of 1 /2 for (e) /Ao, re- 
flecting energy averaging as discussed above (dashed line in 
Fig.|5|l. The cotunneling regularized at resonance (dotted line) 
yields a constant reduction by the small perturbation parame- 
ter Oo: 

8T q + 8T = i3Ao l-4ao 
g^ + gf* 2 1 -205d 

= ^(l-2a )+OK 2 ). (24) 
The full next-to-leading-order calculation (solid line) of S re- 




FIG. 4: Thermopower as a function of j8Ao for Oq = 0.02 at low tem- 
perature. Both the curves for sequential tunneling (gray solid line) 
and sequential plus cotunneling (dotted line) are temperature inde- 
pendent. A full second-order perturbation theory, however, shows 
deviations from this universal scaling behavior due to charging- 
energy gap renormalization. The chosen temperatures are kgT jEc = 
0.01 (dashed line) and 0.0001 (black solid line). We analyze in de- 
tail the reduction of slope at resonance (see inset) and the shift of the 
maximum away from resonance with decreasing temperature. 
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FIG. 5: Average energy of the transported particles as a function of 
temperature for coupling Oq = 0.01. The logarithmic temperature 
dependence typical for Kondo-physics indicates a renormalization of 
the charging-energy gap. We display the sequential-tunneling result 
(dashed line), the "universal" (temperature-independent) result of se- 
quential and cotunneling (dotted) and the full next-to-leading-order 
perturbative result of Eq. <23> (solid). RTA (dot-dashed line) gives 
a similar low-temperature result for the slope in this Figure, but the 
off-set is unknown. 



2. Maximum of thermopower 

The renormalization of the charging-energy gap not only 
modifies the slope of thermopower around n x , it also shifts 
the position of maximum. Figure |6ja) shows the (numer- 
ically determined) position of the maximum as a function 
of temperature. With only sequential and cotunneling taken 
into account, 27 the maximum position approaches a constant 
when lowering the temperature (dotted line) as a consequence 
of the prososed unversal scaling behavior. In a full next-to- 
leading-order theory, however, the maximum position grows 
logarithmically with decreasing temperature (solid line). The 
same low-temperature behavior is reproduced by the resonant- 
tunneling approximation (dot-dashed line). 

Since the exact analytic expression for the maximum posi- 
tion is not transparent, we can gain some more insight in the 
origin of the maximum and the temperature dependence of its 
position by using the following approximate treatment. Close 
to resonance, the average energy from sequential tunneling 
(and thus the thermopower) increases linearly with Aq. Away 
from resonance, the cotunneling contribution dominates, and 
the thermopower decreases proportional to r 2 /Ao. The total 
thermopower including both types of processes is the average 
of the individual thermopower expression, weighted with the 

electrical conductances gy^ cot 



suits in 



-SeT = (e) = 



^Ao/2+^ ot (^r) 2 /A 



scq 

8v ■ 



-8 c v ot 



(26) 



<e>/Ao= 



1 -2oo 2 + 7+ln 



Therefore, the maximum position is roughly at the point 

(25) 

where sequential and cotunneling electrical conductance 
coincide 27 (we disregard any numerical factors here). 



with Euler's constant y = 0.5772 .... The logarithmic tem- 
perature dependence directly reflects the renormalization of 
the charging-energy gap, cf. Eq. (12 1 i . which indicates many- 
channel Kondo physics. 

We also show the result from the resonant-tunneling 
approximation (dot-dashed line), which is obtained 
from numerical integration of Eqs. (|8} and (|9} with 
Eq. H7\ . The logarithmic temperature dependence 
again reflects the renormalization of the level position 
A /A = [1 +2ao(7+ln(j3D/2^))] _1 close to resonance. 
While the logarithmic temperature behavior and consequently 
the slope in Fig.[5]is reliably predicted by RTA, the absolute 
vertical position depends on the choice of the high-energy 
cutoff (here we took D = Eq). 

As thermopower measures the average energy of transport 
it yields a direct extraction of the renormalized charging- 
energy gap A (via the slope of thermopower at resonance). 
This complements in a very appealing manner electrical con- 
ductance measurements, which reveal the renormalization of 
the coupling constant a. As discussed in Ref. 12 and experi- 
mentally observed in Ref s. 1 1 6l and 1 1 71 this renormalization of 
coupling is seen as logarithmic reduction of the maximal lin- 
ear conductance at low temperatures. 




FIG. 6: (a) Position of the maximum of thermopower for ao = 0.01. 
When taking only sequential and cotunneling into account, the max- 
imum's position approaches a constant at low temperature (dotted 
line). In the full next-to-leading-order perturbative calculation (solid 
line) and the RTA (dot-dashed line), the maximum moves away from 
resonance with decreasing temperature, (b) Approximative deter- 
mination of the temperature dependence of the maximum's posi- 
tion [thick solid line as in (a)]. The crossover between sequen- 
tial and cotunneling (squares) becomes temperature independent at 
low-temperature. The correct temperature dependence is reproduced 
when the second-order terms associated with renormalization of the 
system parameters are added to the sequential-tunneling (circles) but 
not to the cotunneling (triangles) contribution. 
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In Fig. |6jb), we show the maximum position determined 
in this way. If sequential and cotunneling processes were the 
only contributions to be considered, the obtained maximum 
position is of the order of the numerically determined value 
and is constant at low temperature (squares). However, in a 
full second-order calculation there are additional terms, as dis- 
cussed above. How do these terms fit into this picture? When 
only looking at the power of the perturbation expansion in Oo 
one might consider these second-order terms belonging to co- 
tunneling. This, however, does result in a completely wrong 
temperature dependence (triangles). The maximum position 
should rather be determined by equating the cotunneling elec- 
trical conductance with that of sequential tunneling plus the 
extra second-order terms, that are interpreted as renormaliza- 
tion corrections to sequential tunneling (circles). In this case, 
the correct temperature dependence is reproduced. This, once 
again, supports the picture of renormalization. 

The main effect of renormalization on the maximum posi- 
tion is that the peak of the electrical conductance gy q around 
resonance is broadened since the renormalized charging- 
energy gap A = A ( )[l — 2ao(const. + \nj3Ec)] is reduced, i.e., 
the system is moved closer to resonance. As a consequence, 
the maximum position moves away from resonance, A max = 
Ao[l +2ao(const. +lnj3£c)] + 0(a ( 2 ), to compensate for this 
renormalization, so that the renormalized maximum position 
is left unchanged. This is indeed the asymptotic behavior 
found in Fig. [6] In conclusion, the temperature dependence 
of the maximum position reflects the temperature-dependent 
renormalization of the charging-energy gap. 

IV. CONCLUSIONS 

In this paper, we investigated the low-temperature proper- 
ties of thermal conductance and thermopower of a metallic 
single-electron transistor. We presented two approximation 
schemes for analyzing higher-order contributions associated 
with quantum fluctuations: A systematic perturbative expan- 
sion in the tunnel-coupling strength and the nonperturbative 



resonant-tunneling approximation. Both these schemes reveal 
qualitatively similar physical effects and coincide for weak 
coupling. In particular, we find that quantum-fluctuation- 
induced renormalizations of the charging-energy gap and the 
tunnel-coupling strength affect the thermoelectric properties. 
They yield logarithmic temperature dependences typical for 
the (many-channel) Kondo effect. 

For the thermal conductance, renormalization of tunnel- 
coupling strength and charging-energy gap results in a sup- 
pression and a broadening of the resonance features, re- 
spectively (see Fig. The effect of charging-energy-gap 
renormalization is most striking in the thermopower. It de- 
stroys the universal low-temperature scaling that would fol- 
low from considering sequential and cotunneling processes 
only (Fig. @}. The reduction of the charging-energy gap is 
reflected in a reduction of the slope of the thermopower at res- 
onance, and in a shift of the maximum's position away from 
resonance (Figs. [5] and |6] respectively). This is due to the 
fact, that thermopower can be interpreted as measure of the 
average energy of transported particles. Therefore, an exper- 
imental observation of the charging-energy gap renormaliza- 
tion in the thermopower would provide an appealing comple- 
ment to the tunnel-coupling renormalization measured in the 
electrical conductance J^ii Certainly measuring such effects 
in thermopower is more challenging, as precise control of 
temperature bias and gate voltage is required. SequentiaU2i2£ 
and cotunneling 21 have been observed in thermopower mea- 
surements some years after these tunneling processes had first 
been investigated in the electrical conductance^^ In contin- 
uing this track record also the renormalization effects on ther- 
mopower predicted in this paper may be within reach of future 
experiments. 
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